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ABSTRACT 

The utility of large Sunyaev Zel'dovich (SZ) surveys for determining cosmological parameters from cluster 
abundances is limited by the theoretical uncertainties in the integrated SZ-flux-to-mass relation, Y-M. We explore 
how non-thermal pressure and the anisotropic shape of the gas distribution of the intracluster medium (ICM) 
impacts Y-M scaling using a suite of hydrodynamical TreePM-SPH simulations of the cosmic web in large periodic 
boxes. We contrast results for models with different treatments of entropy injection and transport, varying radiative 
cooling, star formation and accompanying supernova feedback, cosmic rays, and energetic feedback from active 
galactic nuclei (AGN) and/or starbursts. We find that the gas kinetic-to-thermal pressure ratio from internal bulk 
motions depends on the cluster mass, and increases in the outer-cluster due to enhanced substructure, as does the 
asphericity of the ICM gas (which is substantially more pronounced for the dark matter). The asphericity is less 
dependent on the mass and on variations in the simulated physics. With only a ^ 5 - 10% correction to projected 
(observable) ellipticities, we can infer the 3D ellipticities. We find radii around Rsoq - within which the mean 
density is 500 times the critical density - are the most robust for studying virial properties of clusters, being far 
enough out to avoid the complex "short-distance" physics of the cluster core, having a relatively low non-thermal 
to thermal pressure ratio (^ 20%), and having the smallest variance of gas ellipticity as cluster mass and redshift 
vary. Our simulated K-M-slope roughly follows the self-similar Y ^ M^l^ prediction, except for a steepening 
due to a deficit of gas in lower mass clusters at low redshift in our AGN-feedback simulations. The overall Y-M 
amplitudes with AGN feedback and radiative cooling are lower than for the shock-heating-only case, by ^ 30%. 
AGN feedback enhances slightly the overall F-M-scatter, from ^ 1 1% to ^ 13%, a reflection of accretion history 
variations due to cluster merging. The scatter falls back to ^ 11 % if we select clusters with lower kinetic pressure. 
If we split the cluster system into lower, middle and upper bands of Pkin/Ah, we find a ~ 10% effect on Y-M. 
A 3-split on asymmetry as measured by the long-to-short axis ratio has a < 10% effect on Y-M, but using 3D- 
sphericalized estimates instead of projected (cylindrical) has a a ^ 30% effect. Identifying observable second 
parameters related to internal bulk flows and anisotropy for cluster-selection to minimize F-M-scatter in a (fuzzy) 
"fundamental plane" would allow tighter cosmological parameter constraints. 

Subject headings: Cosmic Microwave Background — Cosmology: Theory — Galaxies: Clusters: General — 
Large-Scale Structure of Universe — Methods: Numerical 



1. INTRODUCTION 

Clusters are the largest gravitationally-collapsed objects in 
the universe, forming at sites of constructive interference of 
long waves in the primordial density fluctuations, the coher - 
ent peak-patches (iBardeen et al.lll986l: [Bond & Myerill996l) . 
The interiors are separated from the Hubble-flow, but main- 
tain contact with the nearby cosmic web through ongoing 
accretion and mergers as they evolve. Although there is a 
strong internal baryon-to-dark-matter density bias, a conse- 
quence of collisional-to-collisionless physics, when cluster- 
scale-averaged the smoothed densities are nearly in the univer- 
sal Hubble-volume-smoothed proportion. Clusters have proven 
to be useful cosmological probes as the rarest collapsed-event 
tracers of the growth of structure in the universe, with a well- 
defined number count that steeply falls as mass and redshift in- 
crease. The number density tail is very sensitive to changes in 
cosmological parameters and primordial non-Gaussianity. In 
clusters most of the baryons are in the form of a hot diffuse 
plasma, the intracluster medium (ICM). The remaining baryons 



are housed in the cluster's numerous stars and galaxies. Obser- 
vations of gas in the cluster system not only reveal the detailed 
astrophysical processes at work in the ICM, but the counts de- 
rived as a function of the global-cluster-observables such as 
thermal energy content can allow for a high precision probe of 
the cosmological parameters defining the count density shape 
and amplitude. 

This scheme of using the cluster system for cosmology must 
rely on simulations capturing the physics at work: the clus- 
ters have been revealed to be too complex for simple spherical- 
ized analytical modelling as the observations have been pro- 
gressively refined and resolutions improved. We hope that the 
basic global observables will be sufficiently robust to the high 
resolution complexities that the cluster/group system can be 
cosmologically useful, but that must be demonstrated by de- 
tailed theoretical work with a necessarily heavy computational 
component. Direct observation of mass or gravitational energy 
or overall binding energy of clusters would be ideal, but we 
are stuck with what can be observed, in the optical, X-ray and 
microwave/radio/sub-mm. Each derived observable from these 
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windows into clusters is fraught with comphcation that requires 
a computational understanding. A thermal Sunyaev-Zel'dovich 
(SZ) (Sunvaev & Zeldovich 1970) probe directly observes the 
integrated Compton-y parameter which is a measure of the 
cluster's global gas heat-energy content, a volume-average of 
the thermal gas pressure, and this is related to gravitational 
energy through the virial relation, so it might be expected to 
provide a robust pr obe. In this paper, and the following se- 
quence. BBPS2 3.4 (lBattaglia etal.ll2011a b c). a follow-on to 
(iBattadia et al.ll20TC we focus on the SZ effect, the Compton 
up-scattering of cosmic microwave background (CMB) pho- 
tons by hot electrons with its unique signature of a spatially- 
varying distortion of the CMB spectrum, a decrement in ther- 
modynamic temperature at frequencies below ^ 220 GHz, and 
an excess above. The SZ signal is proportional to the integrated 
electron pressure, so the hot gas of the ICM dominates the ef- 
fect. The SZ surface brightness is independent of the redshift 
a specific cluster is at. Hence SZ surveys have a different se- 
lection function in redshift and in mass than X-ray and optical 
cluster surveys do, being generically more sensitive to higher 
redshift clusters. The combination of the three probes can pro- 
vide more robust, tighter constraints on cosmological parame- 
ters than any one can alone. 

In a large cluster survey there is a wealth of information 
contained on cosmology and structure formation. The abun- 
dance of clusters, their distribution in redshift, and their spa- 
tial clustering should be determined purely by the geometry 
of the universe, the power spectrum of initial density fluctu- 
ations, and cosmological parameters such as the nns ampli- 
tude of the (linear) density power spectrum on cluster-mass 
scales, erg, the mass-energy density in baryons, dark matter, 
and dark energy, and the equation of state of the latter In SZ 
surveys, the number counts as a function of the total SZ flux 
(integrated Compton-y parameter Y) and redshift and the an- 
gular power s pectrum are two c omplementary p robes of cos- 
mology (e.g., (Birkinshawlll999t ICarlstrom etaI1l2002h . Iden- 
tifying clusters through blind SZ surveys and measuring their 
integrated power spectrum have been long term goals in CMB 
research, and are reaching fruition through, e.^., the South 
Pole Telescope, SPT (e.g ., Lueker et al. 2010; Shi rokoff et al.l 
l20T0t iKeisler et al.l lIoTU [Vanderhndeetal.,,2010(). the Ata- 
cama Cosmology Telescope, ACT (e.g.. Fowler etaH 120101: 
iDunkleyet al l '2010- "Marriage et al." '2010)7 and the Planck 
satellite (e.g., Planck Collaboration et al. 2011b c.a). To de- 
termine cosmological parameters from number counts requires 
understanding the relationship of SZ observables such as total 
SZ flux to fundamental cluster properties such as mass M. And 
to determine them from the SZ power spectrum requires know- 
ing the sum of the squares of pressure profiles of unresolved 
groups and clusters, as well as of resolved ones, weighted by 
the counts. Both depend sensitively on erg, hence can provide 
an independent measure of it. Such extraction from the SZ 
probes is inevitably entangled with the uncertainties in the as- 
trophysical properties of the ICM. This paper deconstructs the 
influence of various physical processes on the Y-M scaling re- 
lation. BBPS2 does the same for the SZ power spectrum. 

Previous work attempted to calibr ate the Y-M scal- 
ing relation through observations (e.g., Benson et al.l l2004t 
Bonamente et al. 2008; Marrone et al. 2009; Andersso n et al.l 
12010; Sayers et al. 201 1; Marrone et al. 201 1), self-calibration 
techniques (e.g., Maiumdar & Mohr 2003, 2004; Lima & Hu 
l2004tlChaudhuri & Ma iumdar 201 1; Nath & Maiumdar 201 1), 
simulations (e.g.. Ida Silva et al. .2004; Motl et al., 2005; 



proac hes (iBodeet al.1 l2007t IShaw et al.l l2008t iMroczkowsO 
1201 Ih . Combining Y-M scaling relations so determined with 
the survey selection function and marginalizing over asso- 
ciated statistical and systematic uncertainties can enable ac- 
curate determination of cosmolo gical parameters. Usi ng a 
small s ample of SZ-clust ers, SPT dVanderlinde et al.l2010h and 
ACT (Seh gal et al.ll201 iT) determined some cosmological con- 
straints, e.g., on erg. However, the errors on erg are dominated 
by systematic uncertainties in the underlying cluster physics, 
making this approach not competitive with other cosmological 
probes. Hence in order to improve upon the determination of 
cosmological parameters, a better understanding of the mass 
proxi es and their scatter is needed (Nagai 2006; Shaw et^ 
l2008t iStanek et al.ll20Tol: I Yang et al.ll2010t iKrause etallBoTTT 

1.1. ICM processes 

Clusters have been increasingly revealed to be complex sys- 
tems as the data has progressively improved, necessitating a 
revision of the simplified pictures popular in the eigh ties for 
interpreting the data. For example, pioneering work bv lKaiseJ 
( 1986) assumed that clusters were self-similar systems with the 
mass determining their ICM thermodynamic properties. As 
shown by subsequent X-ray observations, this self-similar de- 
scription is broken, especially on group scales; low-mass sys- 
tems are l ess luminou s in comparison to the self-similar expec- 
tation (see lVoitll2005l for a review). Studying how non-thermal 
processes such as magnetic fields, cosmic rays, active galac- 
tic nuclei (AGN), star formation, radiative cooling and bulk 
motions contribute to the energy balance and thermodynamic 
stability within clusters is a very active research field. It re- 
mains unclear how these processes vary with cluster radius or 
dynamical state. State-of-the-art simulations are about the only 
tool available for building a consistent picture of clusters. Here 
we contribute to these non-thermal studies by using our simu- 
lations to explore the three effects that influence the Y-M scal- 
ing relation. These are the feedback processes that appear to 
be necessary to explain the thermodynamic characteristics of 
the ICM and avoid a cooling catastrophe leading to too much 
star formation, non-thermal pressure support from bulk mo- 
tions internal to the clusters that are a natural consequence of a 
dynamically evolving structure formation hierarchy, and devi- 
ations from spherical symmetry. 

1.2. Energetic feedback 

In many clusters th e ICM cooling times are much shorte r 
than a Hubble time (iFabianI IT994r ICavagnolo et al.l l2009h . 
which should cause extremely high star formation rates that 
are well beyond what is observed. However, current simula- 
tions with only radiative cooling and star formation excessively 
over-cool cluster centers (e.g ISuginohara & Ostrikeri 119981; 



ISchafer etani2006albn B onaldi et al.ll2007h and analytical ap- 



iLewis et al.l2000tlPearce et al.l2000l) . even with the addition of 
supernova feedback. This leads to too many stars in the clus- 
ter cores, an unphysical rearrangement of the thermal and hy- 
drodynamic structure, and creates problems when comparing 
simulations to observations, in particular for the entropy and 
pressure profiles. Self-regulated, inhomogeneous energy feed- 
back mechanisms by, e.g., AGN are very successful in globally 
stabilizing the group and cluster "atmos pheres", and, in partic - 
ular, preventing the cooling catastrophe ( Churazov et al.l2001h . 
Observations of cool core galaxy clusters show evidence for 
AGN-moderation of the cooling and AGN feedback can po- 
tentially heat the surrounding ICM from kpc sized bubble s 
to hundreds of kpc sized outbursts (iMcNamara et al.l l2005h . 
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In hydrodynamical simulations, it has been shown that incor- 
porating a sub-grid f or AGN fe edback can resolve the over- 
cooling problem (e .g.JSijacki et al. 2007, 2008; Battagli a et alj 
l20T0l:lMcCarflive t al. 2011). The effects of AGN feedback on 
the ICM will mainly alter the cluster and group cores, where 
the actual physics is poorly resolved and understood. These 
effects can s een to be dramatic in X-ray observations (e.g., 
iFabian et al. l 2003), since the emission is proportional to gas 
density squared. Since the SZ signal is proportional to the 
gas pressure, these effects are smaller. Hence, AGN feedback 
should only perturb the integrated thermal SZ signal, with an 
amplitude that is not yet known. 

1.3. Non-thermal pressure support 

Studying non-thermal pressure support from bulk motion in 
clu sters ha s a long history and was first noticed in simulations 
by lEvrardl (119901) who showed that estimates for the bind- 
ing mass of a cluster using a hydrostatic isothermal f3-model 
in comparison to a fit to the surface brightness profile dif- 
fered by 15%. They found that inclusion of velocity disper- 
sion in the hydrostatic isothermal /3-model reconciled this dif- 
ference between binding masses. Including the support from 
residual gas motions in the hydrostatic cluster m ass estima- 
tor im proved the match with the true cluster mass dRasia et aTl 
l2004h. wit h incr easing kinetic pressure at larger cluster radii 
(iLau et al]|2009l) . The amount of energy in these bulk motions 
are of the order of 20% to 30% at radii of in terest for cosmol- 
ogy jBattaglia et a l. 2010; Burns et al.l2010l) . However, kinetic 
pressure support has only recently been included in analytical 
and s emi-analytica l templ ates for the th ermal SZ power spec- 
trum (iShaw et al.l20iaiTrac et al.l201 II) . Of course cosmolog- 
ical hydrodynamical simulations fully include this contribution 
and thus do not require additional modeling of kinetic pressure 
effects. While of importance for correctly interpreting SZ mea- 
surements, the X-ray observations of clusters have been cali- 
brated to partly take this effect into account when determining 
mass from the X- ray inferred tota l thermal energy (e.g., using 
the Yx-M relation. iKravtsov et al .'2006). 

In this paper, our focus is on the effects of bulk motions 
within clusters. These dominate the total kinetic pressure bud- 
get since there is generally a smaller fraction of energy in a 
hydrodynamical turbulent cascade compared to the energy on 
the injection scale which is well resolved for the relevant large 
scale motions. Quantifying turbulence in clusters is becoming 
feasible with sim ulations that include th e modeling of sub-grid 
turbulence (e.g., Ilapichino & Niemeverll2008.) as well as sim- 
ulations wit h magnetohydrody namics and anisotropic thermal 
conduction ( Parrish et al.ll201 II) . 

The method of smoothed particle hydrodynamics (SPH) that 
we are using for solving the inviscid Euler equations in this 
work is perfectly suited for studying large-scale bulk motions 
which dominate the kinetic pressure support due to its La- 
grangian and conservative nature. However, it is known that 
SPH in its standard implementation poorly resolves hydrody- 
namical instabilities, such as of Kelvin-Helmholtz or Rayleigh- 
Taylor type ( Agertz et al. 2007). In non-radiative simulations 
of cluster formation, adaptively-refined mesh codes generate 
a larger core entropy level in comparison to SPH simulations 
which is presumably due to the difference in the amount of 
mixing in SPH and mesh codes and possibly related to a differ- 
ent treatment of vorticity in the simulations (e.g., Frenk et al. 
[19991; IMitchell et al.ll2009t IVazz a et al. 201 1). A recent com- 
parison of a galaxy formation simulation with the SPH tech- 
nique and the recently developed moving mesh code AREPO 



( iSpringel||2010l) enabled - for the first time - to test the na- 
ture of the hydrodynamic solver with otherwise identical im- 
plementations of the gravity solver, the sub-resolution physics, 
and the detailed form of the initial conditions. The moving 
mesh calculations resulted in more disk-like galaxy morpholo- 
gies in comparison to SPH. This difference originated from 
an artificially high heating rate with SPH in the outer parts 
of haloes, caused by viscous dissipation of inherent sonic ve- 
locity noise of neighboring SPH particles, an efficient damp- 
ing of subsonic turbulence injected in the halo infall region 
potentially by the artificial viscosity employed by SPH, and 
because of a higher efficiency of gas stripping in AREPO 
([Vbgelsberger et al. 2011). Based on these results it may be 
questionable to analyse small-scale velocity power spectra and 
2-point structure functions in SPH simulations; however, the 
dominating large-scale bulk motions which is our primary in- 
terest should be followed accurately. This theoretical expec- 
tation is confirmed by analyzing otherwise identical cluster 
simulations run with AREPO and GADGET/SPH that show 
an equal kinetic pressure contribution outside the core region 
(r > 0.05/?20o) for the two numerical techniques (Puchwein & 
Springel, priv. comm., 201 1). 

1.4. ICM shapes 

DM halo shapes have been studied extensively. For SZ ob- 
servations, the shape of the gas distribution of the ICM is im- 
portant, especially in the far field of the intracluster medium 
which contributes sub stantially to the total integrated SZ flux 
(jBattaglia et al.ll20T^ . The assumption of spherical symmetry 
is often made when calculating cluster properties from observa- 
tions and in analytical prescriptions so we would like to assess 
its validity. Semi-analytic models that employ the full three- 
dimensional information of a dissipationless dark-matter-only 
simulation use the shape of the resulting gravitational cluster 
potentials, so it is important to study how such shapes com- 
pare with those that include the dissipational gas component. 
Recent numerical work has shown the impact of cooling and 
star-formation on the propertie s of ICM shape for a sample size 
of 16 clusters (I Lau et af]|201 II) . however, such a study has not 
been extended to a larger sample. Furthermore, the question of 
how energetic feedback in the cluster cores affects ICM shapes 
has not been addressed. 

1.5. Overview 

In this work we explore a large statistical sample of simu- 
lated clusters with identical initial conditions but employing 
different models for sub-grid physics. We quantify the impor- 
tance of non-thermal pressure support and ICM shapes on the 
SZ Y-M scaling relation. In Section |2l we briefly describe the 
simulations and sub-grid physics used. We present our results 
for non-thermal pressure support from bulk motions and ICM 
shapes in Sections [3] and |4] respectively. The impact of these 
processes and changes in the simulated physics on the Y-M 
scaling relation is presented in Section|5] In Section|6]we sum- 
marize our results and conclude. 

2. COSMOLOGICAL SIMULATIONS AND CLUSTER DATA SET 

We d escribed the basic su ite of hydrodynamical models 
used in 'B attagha et alj ( 120101) . We simulate tens of large- 
scale boxes of the cosmic web in order to improve our statis- 
tics of the number of objects while simultaneously aiming for 
a sufficiently high mass resolution to map out the core re- 
gions of those clusters and groups which are the target of cur- 
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rent SZ cluster surveys and which dominate the SZ power 
spectrum signal on scales larger than 1'. Here we charac- 
terize the average behaviour of the properties of the ICM 
over a large mass and redshift range u sing a modified ver- 
sion of the GADGET-2 (I Springelll2005h code which employs 
SPH and treePM for the gravity solver. For each modeled 
physics, we simulate a sequence of 10 boxes of side length 
165 /i"' Mpc = 200 Mpc with periodic boundary conditions, en- 
compassing A^DM = A^gas = 256^ DM and gas particles. This 
gives an initial gas particle mass of Wgas = 3.2 x 10'/i~'Mq 
and a DM particle mass of mum = 1.54 x 10"'/?"' M©. We 
adopt a minimum comoving (Plummer) gravitational smooth- 
ing length of Es = 20 kpc. Our SPH densities are computed 
with 32 neighbours. For our standard calculations, we adopt a 
tilted ACDM cosmology, with total matter density (in units of 
the critical) fi^ = iloM + ^^b = 0.25, baryon density Jib = 0.043, 
cosmological constant Oa = 0.75, a present day Hubble con- 
stant of Hq = lOO/i km s~' Mpc~' with h = 0.7, a spectral index 
of the primordial power-spectrum «s = 0.96 and = 0.8. 

We compare results for three variants of simulated physics: 
(1) the classic non-radiative 'adiabatic' case with only grav- 
itational formation shock heating; (2) an extended radiative 
cooling case with star formation, supernova (SN) energy feed- 
back and cosmic rays (CRs) from str ucture formation shocks 
(for more information on CRs . see iPfrommer et al.l I2006al 
|2S37; EnBHn et al. 2007; JubelgasetalJ|2008|); (3) AGNfeed- 
back in addition to radiative cooling, star formation, and SN 
feedback. Radiative cooling and heating were computed as- 
suming an optically thin gas of primordial composition in 
a time-dependent, spatially uniform ultraviolet background. 
Star formation and supernovae feedback were modelled us- 
in g the hybrid m ultiphase model for the interstellar medium 
of ISpringel & Her nquist (2003). The CR population is mod- 
elled as a relativistic population of protons described by an 
isotropic power-law distribution function in momentum spac e 
with a spectral index of a = 2.3, foUowing lEnBlin et al.l (120071) . 
With those parameters, the CR pressure causes a small reduc- 
tion in the integrated Compton-y parameter (Pfrommeret al] 
120071) . but can result in interesting modifications of the local 
intracluster y-map. 

The AGN feedback prescripti on we adopt for our s tandard 
simulations (for more details see'Battag lia et al.ll2010l) allows 
for lower resolution and hence can be applied to large-scale 
structure simulations. It couples the black hole accretion rate to 
the global star formation ra te (SFR) of the cluster, as suggested 
by [Thompson et al.l ( 120051) . If the SFR is larger than an ob- 
servationally motivated threshold, > 5M0 yr"\ the thermal 
energy is injected into the ICM at a rate which is proportional 
to the SFR within a given spherical region. The AGN feedback 
in these box simulations injects approximately one third of total 
injected energy in the cluster formation phases at z > 2 (analo- 
gous to high-z QSO like feedback), another third in the redshift 
range 1 < z < 2, and the final third below z = 1 (analogous to 
jet/bubble like feedback). These fractions depend moderately 
on the the numerical resolution; increasing the resolution en- 
ables to resolve the growth of smaller halos at earlier times and 
cau ses a higher fraction o f energy injection at higher redshifts 
(see lBattagha et alj|2010l for a discussion). 

We define the virial radius of a cluster, Rj^, as the radius 
at which the mean interior density equals A times the criti- 
cal density, Paiz) (e.g., for A = 200 or 500). For comparison, 
we will use an alternative definition of the virial radius, /?A.m, 
where the mean interior density is compared to the mean mat- 



ter density, Pmiz). For clarity the critical density and the mean 
matter density are, 

P.Az)=^[^u.il+Z? + ^A\, (1) 

Pu.{z)=^^M+zf. (2) 

Here we have assumed a flat universe (i^m + fiA = 1 ) and are 
only interested in times after the matter-radiation equality, i.e., 
the radiation term with 57i is negligible. We chose to define the 
virial radius with respect to the critical density in continuity 
with recent cluster measurements. The merits and utilities of 
both these definitions are discussed later in AppendixICl 

We apply the following two-step algorithm to compute the 
virial mass of a cluster in our simulations. First, we find all 
clusters i n a given sna pshot using a friends-of-friends (EOF) al- 
gorithm (iHuchra & G eller 1982). Then, using a spherical over- 
density method with the EOF values as starting estimates, we 
recursively calculate the center of mass, the virial radius, R/\, 
and mass. Ma, contained within R/\, and compute the radially 
averaged profiles of a given quantity with radii scaled by /?a- 
We then form a weighted average of these profiles for the entire 
sample of clusters at a given redshift unless stated otherwise. 
We use the integrated Compton y-parameter as our weighting 
function, 

Ya = ^ I ^ P,(r)ATTr'dr « E,^{< R^) , (3) 

where is the Thompson cross-section, me is the elec- 
tron mass and is electron pressure. For a fully ionized 
medium of primordial abundance, the thermal pressure P = 
fe(5XH + 3)/2(XH + l) = 1.932Pe, where = 0.76 is the pri- 
mordial hydrogen mass fraction. 

3. NON-THERMAL CLUSTER PROFILES 

Several simulations dEvrardl 119901; iRasia et alj 120041; 
|Lau et al. 2009) showed that the kinetic pressure from bulk 
motions contributes a small but still significant amount of en- 
ergy within 7^500 and this importance increases for larger cluste r 
radii dLau et al.ll2009t iBattaglia et alllloTot iBurns et alJiMoh . 
Thus, it is important to accurately quantify the kinetic pres- 
sure contribution as it biases the hydrostatic cluster masses and 
is significant to the total energy budget within clusters. There 
are two kinetic pressure contributions, namely large-scale, un- 
virialized bulk motions and subsonic turbulence. For a Kol- 
mogorov power spectrum of turbulence, the energy is domi- 
nated by the largest scales which we resolved and characterize 
in our simulations. Hence we believe that our approach cap- 
tures the majority of the kinetic pressure contribution. 

3.1. Kinetic pressure support 

The internal bulk motions in the medium can be quantified by 
the mass-averaged velocity fluctuation tensor, {6ViSVj), which 
is associate with the kinetic pressure (stress) tensor, 

Pi^n,ij=p{SViSVj), 

PMn = TTP^j3 = p{dV-SV)/3, 

5V = a(v-v) + aH(z)(x-x) . (4) 

In this paper we focus on the trace, i\in, which we refer to 
as the kinetic pressure. Issues associated with the anisotropic 
stress tensor will be made explicit where they appear, and are 



ON THE CLUSTER PHYSICS OF SUNYAEV ZEL'DOVICH SURVEYS I 



5 




Figure 1. The ratio of kinetic and tliermal pressure support, Pkm/Pth, depends on mass and redshift. We show the median of /\i„/P|h as a function of radius for 
the AGN feedback simulations for various mass bins at z = (left) and as a function of redshift for a fixed mass bin (right). We additionally show the 25''' and 75"' 
percentile values for the lowest mass bin aX z = (dotted). In both panels we illustrate the 1 and 2 a contributions to Y/\ centered on the median for the feedback 
simulation by horizontal puiple and pink error bars which extends out to 4S2oo iBattaglia et al. 2010). Two analytical models for the Pi^i,, by Shaw et al. (2010) and 
ITrac et atl <201 ll) are shown with the dash dot and dashed lines, respectively. The iShaw et aU (12010 ) model matches our result in the mass bin 2.7 X IO'^'Mq < 
M200 < 4.2 X 10''* M0 at intermediate cluster radii (this mass bin best represents the mean mass of their sample at redshift zero), but also illustrates the need for a 
mass dependence in future analytical models. The dependence of Py„/P|i, on cluster mass is driven by the variation of Pi;;,, with mass (see Fig. |2]below). 




Figure 2. Mass and redshift dependence of Pth and P^cm, normalized to an empirical fit from BBPS2, Pflt, to the scaled thermal pressure, Pth/PA- We show the mean 
scaled thermal and kinetic pressure profiles at z = as a function of radius for the AGN feedback simulations in vaiious mass bins (left), and for vaiious redshifts at 
fixed mass bin (right). 



also explored in more detail in BBPS3. The code uses comov- 
ing peculiar velocities, which are translated into the internal 
cluster velocities relative to the overall mean cluster velocity in 
the Hubble flow by the relation given, where H(z) is the Hub- 
ble function, a is the scale factor, v (= dx / dt) is the peculiar 
velocity and x is the comoving position of each particle. The 
gas-particle-averaged cluster bulk flow within /?20() is "0 and the 
center of mass within /?2()o is x. 

The radial profiles of the kinetic-to-thermal pressure, 
Pkin/Pth, shown in Fig. [1] for various mass bins demonstrate 
an overall mass dependence at all cluster radii, predominately 
driven by the variation of P^n and not Pth with mass. We show 
this explicitly in Fig. |2]where we scale Pjdn with the virial ana- 



logue of the thermal pressure. 

Pa = GMAApcr(z)/b/(2PA), /b = ^b/^,^- (5) 

This behavior reflects the average formation history of galaxy 
groups and clusters which, according to the hierarchical pic- 
ture of structure formation, sit atop the mass hierarchy, with 
the most massive clusters forming and virializing near to the 
present time. In contrast, the median galaxy group (M200 = 
IO'-'Mq) has stopped forming today as can be seen by the 
dramatically decreasing mass accretion rates implying that the 
associated virializing shocks have dissipated the energy asso- 
ciated with the growth of th ese objects and hence decreasing 
the kinetic pressure support dWechsler et alj 120021 : IZhao et alj 
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I2OO9I: iPfrommeretan IMlh . The semi-analytic model for 

falls 



analyt 

non-thermal pressure support by Shaw et alj (1201 Oh talis m 
the middle of the mass bins chosen since this model results 



from a sample of 16 high resolutions adaptive mesh r efinement 
(AMR ) simulations of individual galaxy clusters (iLau et alJ 
l2009h which have a similar mass range. We provide a simple 
fit for the mass dependence of Pkin/^'th in Appendix|A] We find 
that the radius at which Pyn = Ah is just beyond the spherical 
collapse definition for 7?vir from lBrvan & NormanI (1 19981) . 



Rvh = 



3Mv 



47rA„(z)pcr(z) 



1/3 



(6) 



where A„(z) = 18tt^ + 82[i^(z)- l]-39[rj(z)- l]^ and n(z) = 

^m(i + z)^ [flmi^ + z)^ + flA] ■ Hcncc, this radius represents a 
possible physical definition for the virialized boundary of clus- 
ters. 

The redshift evolution of Pkin/Ah is dramatic. At higher red- 
shift, fkin is increasing faster than Pth over all radii (cf. Fig. 
|2|, such that at z = 1, Pkin/^th is approximately twice that at 
z = 0. In the picture of hierarchical structure formation, at any 
given redshift the most massive objects are currently assem- 
bled and hence show the largest kinetic pressure contribution in 
comparison to smaller objects that formed on average earlier. 
Or equivalently, at fixed cluster mass, the relative contribution 
from kinetic pressure and the relative amount of substructure 
increases with redshift. In particular, the relative mass accre- 
tion rates increase from z = 0toz = 2bya factor 3 for clus- 
ters (M200 = l O'^Mn) and 10 for groups (M200 = lO'^M©) (see 
iPfrommer et al. 2011; Gottloberet al. 2001). This strong evo- 
lution in Pkin/Pth is lessened by a different choice of scaling 
radius, i.e., if we normalize by ^200,m instead of R200 (cf. Ap- 
pendixlcl- Although this ratio cannot be observed, we will use 
it as an indicator for the dynamical state of clusters in our simu- 
lations. Results from Lau et al. (2009) find a similar correlation 
between Pyn and the X-ray definition of dynamical state, frorn 
a smaller sample of 16 clusters. At z = 1, the Shaw et al J (1201 Oh 
semi-analytic model for non-thermal pressure support does not 
match our simulations as well as it does at redshift zero0 

The formation of clusters and the associated accretion of sub- 
structure are driven by the depth of the cluster gravitational po- 
tential. Therefore, it is not surprising that we find kinetic pres- 
sure support to be ubiquitous in the three differently simulated 
physics cases (cf. Fig. [3]l. Looking at the median of this non- 
thermal pressure support we find similar radial profiles within 
the 25* and 75* percentiles of the complete distribution of 
clusters. In the AGN feedback simulations we find marginally 
lower values for Pkin/Ah- These differences are well within 
the 25* and 75* percentiles implying consistency across dif- 
ferently modeled physics. Thus, our model of AGN feedback 
does not significantly alter the kinetic pressure support at low 
redshift although there seems to be a hint that this may be the 
case at larger radii at redshifts z ^ 1 which approach the peak 
of the AGN luminosity density. 

We have so far focused on the trace of the kinetic pressure 

' Our kinetic pressure contribut ion is larger at the center compared to that 
in the model by Sha w et alj 1201(1 ). This discrepancy is probably a manifes- 
tation of the well-known core entropy problem in numerical simulations. In 
(adaptive) grid codes there is a larger level of core entropy generated in com- 
parison to SPH codes implying that the enhanced entropy (which results from 
dissipating gas motions) is accompanied by a smaller amount of kinetic pres- 
sure. This is presumably due to the difference in the amount of mixing in SPH 
and mesh codes and possibly related to a different treatment of vorticity in the 
simulations (e.g.. lFrenk et alJ1999l : lMitchell et al.l2009tlVazza et alJ2011h . 



1.7x 10'*M®<M,„<2.7 X lO'^Ms 

AGN feedback 

Radiative cooling 

Sliock iieating 

; = 




Figure 3. Shown is the median of Pkin/^th as a function of radius for differ- 
ent physics models at z = (solid) and z = 1 (dashed) with the 25"" and 75"" 
percentile values shown for the AGN feedback simulations at z = (dotted). 
Results are shown for the mass bin 1.7x IO'-^Mq < Mjoo < 2.7 x lO" M© to 
take out the dependence on mass of Pkin/^th- The kinetic pressure contribution 
is similar for our differently simulated physics, suggesting that gravitational 
processes dictate that contribution (while AGN feedback slightly decreases the 
kinetic pressure contribution, especially for higher redshifts). The horizontal 
purple and pink error bars have the same meaning as in Fig.[T] 
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AGN feedback, z = 

1.1 X 10'*Mg,<M2„< 1.7x lO^Mg 

1.7x 10"'Mg,<Mj„<2.7x 10'* Mg 
2.7 X lO" M5,< Mj„ < 4.2 x lO'"* Mg 

4.2 X lO" M®< A/j„ < 6.5 X 10" 

6.5 X lO" Ms< A/j„„ < 1.01 X 10" Ms 

1.01 X 10" Mo< */,„„< 1.57 X 10" Ma 




Figure 4. Assessing the bias in hydrostatic masses, Mhse> due to the kinetic 
pressure. The median of Mhse/Moi as a function of radius for AGN feedback 
simulations for various mass bins, with the 25"' and 75"" percentile values 
shown for the smallest mass bin (dotted). Assuming hydrostatic equilibrium 
for all clusters of a given mass will bias the mass values low by 20 to 25%. 
The scatter about the median - represented by the 25"" and 75"" percentiles 
- amounts to approximately 5%. This bias is not representative for a relaxed 
cluster sample which will likely have a smaller bias since th e calibration of 
such a sample against numerical cluster simulations shows IKravtsov et alj 

tensor. We treat in detail the velocity anisotropy in BBPS3. 
The main results are that the core regions near the center, the 
velocity distribution starts to become isotropic for the gas in 
groups and (to a lesser extent) for the DM and gas in larger 
clusters. The positive values of velocity anisotropy around the 
virial radius indicate (radial) infall, whereas the strong decrease 
at even larger radii (very noticeably in the DM) is caused by the 
turn-around of earlier collapsed shells, which minimizes the 
radial velocity component such that the tangential components 
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dominate the velocity. 

3.2. Hydrostatic Masses 

Even if clusters are in hydrostatic equilibrium, balancing the 
gravitational force to the pressure gradient yields 

VP = pg ^ -pGM(< r)r/r^ for spherical symmetry, (7) 

and using this relation to estimate cluster masses will give the 
wrong result if one does not include non-thermal pressure, in 
particular the kinetic pressure. 

As other s have shown (Evrar3 Il990t iRasia et alJ 12004 
[Lau et al] 120091) . assuming that all the pressure in Eq. (7} is 
thermal (P = fth) is incorrect; for clarity we define Mhse to be 
the mass derived using P = P^h- Comparing Mhse to the true 
mass inside a given radius, Mtot, we find that Mhse on aver- 
age underestimates Mtot by 20-25% depending on the radius 
(cf. FiglH. This bias is almost independent of cluster mass out 
to /?5oo, the current maximum radius typically observed by the 
X-ray telescopes Chandra and XMM-Newton. We can under- 
stand this weak mass dependence by rewriting the total pres- 
sure P = P{P,b/P+Pkin/P)- Since Pkin/P cx Mj^q (see Fig. [HI, 
the hydrostatic mass estimates inherit a similarly weak depen- 
dence on mass. Thus, an overall correction to the hydrostatic 
mass is reasonable for these measurements. 

Individual clusters can stray from this generalization, since 
each cluster has a unique dynamical state and formation his- 
tory. These deviations are suggested by the scatter of ^ 5% be- 
tween the 25* and 75* percentiles of the complete distribution. 
For cluster samples that are selected against major mergers (for 
which the assumption of spherical symmetry will also be ques- 
tionable), the correction factor will necessarily be smaller, e.g., 
for quality X-ray data of a Chandra sample, the hydrostatic 
mass correction was f ound to be of the order Mhse ~ 10-15% 
(iKravtsov et al.l[2006l) . 

4. CLUSTER SHAPES 

Generally, we expect clusters to be triaxial since they grow 
by accretion and through merging along filamentary structures 
that impose tidal gravitational forces upon the forming clusters. 
Following Dubinski & Car lbergI (1199 1). we estimate this non- 
sphericity of cluster gas and dark matter (DM) by computing 
the normalized moment-of-inertia tensor, 

{5x5x).^Xr<R\w)=^ (8) 

for several weightings Wq of the contribution of particle a that 
lies within a given radius R. The tensor measures the variance 
in the spatial fluctuations within R, with 5xi = the devia- 
tion of the particle position x, from the centre-of-mass x,(< R) 
of the region. Using mass weighting, Wa = for the dark 
matter or the SPH gas particles gives the moment-of-inertia in 
its usual form. It has an effective p{x)x^d\nx reach in its probe 
of the unit vector combination 5xi6xj, hence preferentially feels 
the outskirts, near R. In the Appendix |B] we explore how our 
results are modified with a weight = Aq less sensitive 
to the outskirts: this weight just mass-averages the unit vector 
product, hence emphasizes the more isotropic interior. Since 
the tSZ signal is of primary interest to us, we also consider 
thermal-energy weighting, with Wq, = m^Ta the product of the 
mass and temperature. 

We quantify the asphericity of cluster gas and DM by two 
parametrizations: the axis ratios (in particular the ratio of the 



largest-to-smallest main axis, c/a); and the three-dimensional 
asymmetry par ameters for symmetric tensors introduced by 
BBKS (Bardee n et al.l ll986l). Both use the eigenvalues A, of 
the moment-of-inertia tensor at a prescribed radius, ordered by 
Al < A2 < A3. The ellipsoid associated with the tensor has axis 
length s a = \/A7, b = -y/Aa, and c = \/X; ■ (We note that lLau et alj 
(I2OI Ih express their results in terms of eigenvalue ratios, defin- 
ing a' = Xi, b' = X2, and c' = A3). The eigenvectors £■,■ associated 
with A, are also used to rotate the clusters to their principal axes 
and to explore alignment variations with radius. 
The BBKS-style asymmetry parameters are defined by 

A1-A3 

A] — 2A2 + A3 
P = ^ ' 

A = Ai + A2 + A3. (9) 

Following BBKS (iBardeen et al.ll 19861) . we refer to e as the el- 
lipticity and p as the prolaticity. When p is positive the clusters 
are prolate, and when p is negative the clusters are oblate. Thus 
we also define the oblateness o = -p6(-p), equal to \p \ for neg- 
ative p and zero for positive p. Although e is the most strik- 
ing indicator of cluster elongation, the degree of prolateness or 
oblateness are also necessary to specify the general triaxiality 
of the morphological configuration on a given smoothing scale. 

4.1. Overall shapes and their profiles 

We rotate all clusters into the moment-of-inertia tensor frame 
using the eigenvector matrix E, so x' = Ex. The output order- 
ing is arbitrary; we choose the convention that the major axis is 
aligned with the x-axis and the minor axis is aligned with the z- 
axis. In Figs.|5]and|6]we show the results for DM and gas which 
have been obtained by computing the weighted moment-of- 
inertia tensor within 3/?200, rotating into the moment-of-inertia 
tensor frame, and stacking the respective distributions, i.e., gas 
density and pressure as well as DM density. The rotated con- 
tours show obvious elongations along the major axis; with the 
ellipticity being larger at z = 1 in comparison to z = 0. The elon- 
gation is larger for the DM distribution in comparison to the 
gas density and pressure which show very similar behaviour. 
Even in the rotated stacked distributions, the innermost contour 
lines become more spherical because they are intrinsically less 
elliptical (see below) and because the main axes of the inner 
distributions are twisted relative to those at 3/?200 so that their 
ellipticity partially averages out to become more spherical (see 
Sect.O. 

In order to quantify these results, we show the mass depen- 
dence and redshift evolution of the ellipticity within /?2()() in 
Fig. |7] Due to the dissipationless nature of DM, its elliptic- 
ity is larger (smaller ratio of c/a) in comparison to that of 
the gas. This is because the kinetic energy of the accreted 
gas is dissipated at cluster accretion shocks - a process that 
erases part of the memory of the geometry of the surrounding 
large scale structures and their tidal force field. Those accre- 
tion shocks are typically forming at radii > /?200 as suggested 
by numerical simulatio ns (Miniati et al. 2000; Rvu et al. 20031; 
iPfrommer et al.l2006bl;[SkiUman et al..2008;, Vazza et aL2009l) 
or indirectly by the action of shock waves on radio plasma 
bubbles, which represents a novel method of finding forma- 
tion shocks (e.g., Ensslin et al. 2001; Pfrommer & Jones 20ll]). 
Following these qualitative considerations, it is not surprising 
that the ellipticity of the gas distribution does not show any 
mass dependence while the DM distribution of more massive 
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Figure 5. Stacked density and pressure distributions with and without rotations into the principle axis frame of the correspondingly weighted moment-of-inertia 
tensors at z = 0. Left: We compare rotated distributions of the gas density (red) and pressure (blue) to the non-rotated stacked gas density (black) at z = 0. Right: Shown 
is the same as on the left for DM. The non-rotated clusters average out to form spherical iso-density contours, while the rotated clusters clearly show elongations 
along the major axis (defined here as the x-axis). The thicker lines approximately show the radii i?2500, ^500 and R200 from the inside out. These contours have been 
smoothed to a pixel size of 0.09S200- The horizontal purple and pink error bars have the same meaning as in Fig. [T] 




Figure 6. Same as in Fig. \5\ but at z = 1 . 
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Figure 7. We show axis ratios of clusters that are obtained by computing the moment-of-inertia tensor of the gas (red) and DM mass distributions (blue) within R200 
and stacking those in bins of cluster mass and redshift. The resulting mean and standard deviation of the axis ratio c/a is shown as a function of M200 at z = (left 
panel) and at fixed average mass bin of M200 = 2.4 X 10*** /;'' Mp) as a func tion z (right panel). See Table[T]for fit values; here, we have chosen to quote Mq 
to compare directly with the dissipationless simulations bv lKasun & EvrarH 12005. KE05). Shocks dissipate the kinetic energy of the gas which causes larger axis 
ratios/smaller ellipticities, this couples through gravity to the DM distribution and sphericalize their axis ratios, resulting in smaller ellipticities in comparison to 
dissipationless simulations alone, e.g., by KE05 (that do not follow the hydrodynamics of the gas). 




Figure 8. Average cluster axis ratios and ellipticities for the DM mass (blue), gas mass (red) and pressure (green) distributions within a scaled radius r/RjoQ. Left: 
Shown is the axis ratio c/a as a function of scaled radius for all simulated physics models. In the bottom panel we show the relative differences of the shock 
heating (long-dashed) and the radiative cooling simulations (short-dashed) with respect to AGN feedback simulations. Right: We show the ellipticity, prolaticity and 
oblaticity as a function of scaled radius. The bottom panel shows the relative differences between the pressure and gas density weightings of the moment-of-inertia 
tensor. The axis ratio c/a and ellipticity show the same trends. We find clusters to be more prolate then oblate. In the regions beyond Ryj, the sudden decrease in 
the a xis ratios can be attributed to other nearby groups and collapsed objects (that is also seen as an enhanced density clumping at these radii, INagai & LaiJI20lTI : 
IBatta glia et al. 20 lie). The pressure- weighted shapes tightly track the density-weighted shapes with deviations of less than 5%. The horizontal purple and pink error 
bars have the same meaning as in Fig. [T] 
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Table 1 

Axis ratio fits for cluster as a function of mass and redshift. 





Bm 


Am 


Bz 


Az 


DM 


0.665 ±0.009 


-0.031 ±0.006 


0.692 ± 0.009 


-0.12 ±0.04 


KE05 


0.631 ±0.001 


-0.023 ±0.002 


0.652" ±0.001 


-0.086 ± 0.004 


Gas 


0.87 ±0.02 


-0.003 ±0.01 


0.88 ±0.03 


-0.08 ±0.12 



" We use a re-normalized value from M200 = 1x10/? Mq 
toM2oo=2.4x lO''*r'M0. 



clusters shows a larger ratio of c/a in comparison to smaller 
systems. However, the ellipticity of the gas and DM distribu- 
tion are increasing as a function of redshift, at about the same 
rate. This can be understood by the fact that 1) a given mass 
range of clusters shows a larger degree of morphological dis- 
turbances/merging at higher redshifts which probe on average 
dynamically younger objects and 2) the redshift evolution of 
the velocity anisotropy (cf. BBPS3) which shows that the av- 
erage location of accretion shocks moves to smaller radii (if 
scaled by /?2{)o)- Hence at larger redshifts, also the gas distribu- 
tion probes the infall/pre-accretion shock region that is shaped 
by the tides exerted by the far-field of clusters. 

We co mpare the results from o ur simulations directly with 
those of iKasun & Evrard' ('2005) in Figure |7] and Table [U 
Other work (Allgoo d et al . 2006; Gottlober & Yepes 2007; 
iMaccio et all UOO^ ILau et al l [20111) on DM and gas shapes 
have used different mass definitions, axis definitions and cos- 
mologies than we do, making quantitative comparisons diffi- 
cult, but we note that these various results were shown to be 
consistent with Kasun & Evrard (2005). For the mass and red- 
shift functional fits, Kasun & EvrardI (120051) define c/a(M) = 
Bm(1+Am ln[M/10i5/i-'Mo]) and c/a(z) = B^(l+z)^\ respec- 
tively. The axis ratios that we find for the DM mass de- 
pendence and redshift evolution have slopes consistent with 
iKasun & Evrard (2005), but as an overall trend, our axis ratios 
are more spherical than theirs. Some of this may be traced to 
differing cosmological parameters, but some may be because 
the less aspherical baryons may have an impact on the DM el- 
lipti city. The ef fect of baryons on DM has been explored in, 
e.g.. lRudd et"al] (l2008). 

In the following, we will show radial profiles of the axis ra- 
tios and asymmetric parameters that are obtained by comput- 
ing the moment-of-inertia tensor at 30 different radii for each 
cluster. In Figure [8] we report on the overall radial distribu- 
tion of c/a and ellipticity in the gas and DM distributions. 
Within R200, the ellipticities of gas density and pressure are 
rather flat at a level of c/a ~ 0.85-0.9. As laid out above, 
this is because dissipation effects at the accretion shocks cause 
an effective sphericalization and erase the memory of large- 
scale tidal fields. In contrast, ellipticities are increasing for the 
DM as a function of radius due to the dissipationless nature of 
DM, i.e. c/a decreases from values around 0.8 in the center 
to 0.7 at /?20()- The radial behaviour may be due to increased 
tidal effects on DM substructures at small radii which causes 
a dram at ic drop of their ce ntral mass density dSpringel et alj 
l2008allbl:|Pin"zke et al.ll20TTb . Effectively this causes a redistri- 
bution of a clumped (elliptical) to a smooth distribution that 
is able to couple more efficiently to the (more spherical) gas 
distribution. Studying the asymmetric parameters, we find that 
if a cluster is prolate, it is on average more elliptical than an 
oblate one that is always close to spherically symmetric. 

We find that the average axis ratios and ellipticities have a 
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Figure 9. Top: We compare the average of the 2D axis ratios (dashed) of three 
random orthogonal projections to the 3D axis ratios (solid) for the DM mass 
(blue), gas mass (red) and pressure (green) from the clusters in the AGN feed- 
back simulations. Additionally, we show the linear correlation coefficient, rs 
(dotted fine), between the projected 2D and the 3D axis ratios. The horizon- 
tal purple and pink error bars have the same meaning as in Fig. [T] Bottom: 
Shown is the relative differences between the projected 2D and 3D axis ratios. 
While the relative difference between the projected 2D and 3D axis ratios vary 
between 15-20% for the DM mass distribution (with the 2D axis ratios being 
more spherical), the relative differences are smaller for the mass and pressure 
distribution of the gas, with values between 5-10%. As expected, the projected 
2D and the 3D axis ratios are correlated with an increasing correlation coeffi- 
cient at larger radii which suggests that the substructure distribution that drives 
the asphericity also causes this correlation. 

pronounced break in their slopes at r ^ 1.5/?200- The break in 
the ellipticity arises from substructure in the cluster outskirts. 
Recent X-ray observations of the Perseus cluster find a strong 
signature of clumping in gas density (Simionescu et al. 201]]); 
qualitatively consiste nt with the fin dings in simulations but not 
quantitatively (Nagai & Laull20in) . This gas density clump- 
ing is a direct tracer of substructure and becomes important at 
roughly the same radius where we find the break in the elliptic- 
ity. Interestingly, this effect is not only seen in DM and gas but 
also in pressure, which suggests that the pressure is clumped 
in a similar fashion as the gas density. In order to accurately 
model the outskirts of clusters, semi-analytic models will need 
to properly deal with the substructure. In the Appendix iBl we 
show that one can attempt to counteract or lessen the impact of 
substructure on the gas, pressure and DM shapes by including 
an r"^-weighting when calculating the moment-of-inertia ten- 
sor (cf. Eq. (O). In future work, we will further explore the 
issue of substructure. 

4.2. How shape profiles depend on modeled physics 

We also address the influence that changes in the simulated 
physics has on cluster shapes in Fig. [8] While the ellipticity 
of the gas is slightly larger in non-radiative models, it is very 
similar for the gas distribution in our radiative models (radia- 
tive cooling and star formation with and without AGN feed- 
back). Dissipating accretion shocks seem to explain the over- 
all behavior rather well and the different physical models only 



ON THE CLUSTER PHYSICS OF SUNYAEV ZEL'DOVICH SURVEYS I 



11 



marginally change the cluster shapes in the gas. In the DM, 
however, there is still a pronounced difference among our two 
radiative physics models with the ellipticities of the AGN feed- 
back model being larger that in our pure radiative model. This 
small ellipticity is a remnant of overcooling that our pure radia- 
tive model suffers with an associated star formation rate that is 
unphysically high. Most of these stars form out of the cold, 
dense gas in the core region which causes a decreasing central 
pressure support so that gas at larger radii moves in adiabat- 
ically and causes a deeper potential which in turn causes the 
DM to adiabatically contract. Enhanced dissipation processes 
in the gas sphericalize the potential which is then communi- 
cated to the DM during this central settling. We find that in- 
cluding AGN feedback counteracts the overcooling issue and 
modifies the DM shapes on the level of 5% in comparison to 
our pure radiative simulations (cf. FigO. 

Our general trends are similar to those reported by lLau et al.l 
(1201 Ih who also find that the DM distribution is more spher- 
ical for radiative simulations in comparison to non-radiative 
models. However, the differences between radiative and non- 
radiative simulations are not as extreme as those found in 
iLauet all (12011.) . since our radiative simulations do not have 
the level of (catastrophic) cooling in the central regions, be- 
cause their simulations have higher resolution and include 
cooling from metals which we do not. The AGN feedback 
stabilizes the cooling and, thereby, lessens this sphericalizing 
effect on the DM ellipticity. 

4.3. Projected and intrinsic shapes 

In order to tie the underlying 3-dimensional structure of clus- 
ters to observable 2-dimensional projections, we compare the 
intrinsic 3D axis ratios to axis ratios of random 2D projections, 
i.e. we project the DM density and gas density/pressure distri- 
butions along a randomly chosen direction and then compute 
the 2D moment-of-inertia tensor. The results are shown in Fig- 
ure |9] We find that the 2D axis ratios {b/a) for both, the gas 
density and pressure are systematically closer to unity than the 
3D ratios c/a by ^ 5-10% (for the virial region and the cen- 
tral part). In the case of the DM distribution, the projected 
(2D) axis ratios are on average a ^ 15% underestimate of the 
intrinsic (3D) axis ratios. Using the linear correlation coeffi- 
cient statistic (rs), we find that the random 2D axis ratios are 
strongly correlated with the intrinsic 3D axis ratios with a ra- 
dially increasing correlation strength (cf. Fig|9]l. As expected, 
c/a serves as a limit to the observable projected (2D) axis ratio. 
We find that the mean 2D axis ratio, as a function of cluster ra- 
dius, closely tracks c/a, modulo a roughly constant ~ 5 - 10% 
bias (for the gas density and pressure). 

4.4. Mass and redshift dependence of shape profiles 

Both the density- and pressure-weighted ellipticities show 
the same general trends with radius and cluster mass. The ellip- 
ticity increases with increasing cluster mass by ^ 50% over the 
mass ranges shown (cf. FigfTOli. On the right-hand side of Fig. 
[Tol we show the redshift evolution of the cluster shapes and 
find that the ellipticity is a stronger function of redshift than the 
mass. For increasing redshift, the break in the ellipticity profile 
moves to smaller radii (when scaled to 7?2oo)- Both behaviors, 
the mass and redshift dependence can be understood in the hi- 
erarchical picture for structure formation, where clusters show 
increased mass accretion rates and hence an increased level of 
substructure for larger clusters (at a given redshift) or, equiv- 
alently, for a cluster of given mass at higher redshifts which 



probe on average systematically younger systems. Similar to 
the non-thermal pressure support, the redshift evolution found 
in the ellipticities are lessened by a different choice of scaling 
radius (cf. Appendix ICli. This result suggests that using a sin- 
gle (constant) ellipticity profile for clusters is not sufficient for 
percent level accuracy. 

Pressure-weighted ellipticities are marginally more spherical 
than the density-weighted ellipticities for r < R^m (cf. Figs. [8] 
and [Toll. However, between TJgoo and 27?g()o the behavior is re- 
versed. This is because the core region shows a smaller kinetic 
pressure support implying that hydrostatic forces had time to 
act and to smooth out the pressure distribution whereas at larger 
radii, pressure-weighted ellipticities are affected more by in- 
fall caused by a noticeable pressure clumping at these radii (cf. 
BBPS4). We have shown that scales around /?5oo are the most 
robust for studying the virial properties of clusters, since all 
ellipticities show only ~ 10% redshift evolution and mass de- 
pendence on these scales. Taken together with the relatively 
modest degree of non-thermal pressure support on these scales 
found in §Q]and that these scales are far enough out to avoid the 
complications of the intricate "short-distance" physics of the 
cluster core region, we can give further justification for what 
is already the practice in the X-ray cluster community, driven 
by the nature of the X-ray data, namely a focus on /?5oo and 
environs for "global" cluster properties. 

4.5. Alignment variations and semi-analytical models 

Semi-analytic models for the baryon distribution in clusters 
include an underlying assumption that baryons will arrange 
themselves along equipotential surfaces (or in some cases the 
DM density-weighted surfaces). Given the importance of this 
assumption, we test its validity in our simulations. In Figure]?] 
we plot the ratio c/a for both dark matter and gas as function of 
cluster mass and redshift. W hile c/a for DM ha loes decreases 
with halo mass as expected (*Jing & Sutoll2002h . we find that 
c/a is constant for the gas distribution. This i s potentially a 
problem for semi-analytic models of ICM gas (' Ostriker et al.l 
2005; Bode et al. 2009), which solve for the resulting gas dis- 
tribution in a DM-dominated gravitational potential as obtained 
from dissipationless simulations while allowing for a constant 
non-thermal pressure contribution (in the latter case). How- 
ever, the gravitational potential from the DM is more spherical 
than the underlying matter distribution (e.g. Lau et al. 20111), 
and so the semi-analytic shape estimates are not as discrepant 
as one might expect from Figure |7] 

A more important issue is the alignment of the gas or pres- 
sure with respect to the DM. We calculate the angular differ- 
ence between the major axes of the DM and those of the gas 
and pressure major axes at a given radius, using the moment- 
of-inertia tensor eigenvectors £'i_dm('") •■E'i,gas('")- When calcu- 
lating misalignment, the major axes of nearly spherical objects 
are poorly defined quantities. To avoid this problem we cal- 
culate, in each radial bin, a weighted average using 1 - c/a as 
the weight. Furthermore, we exclude the region inside 0.37? 200 
since the gas and pressure shapes are nearly spherical, with 
c/fl|DM > 0.75. On average at a given radius, the cluster gas 
and pressure are 20-30 degrees misaligned from the major axis 
of the DM (cf. Fig. fTTT i. In the next section we show SZ mea- 
surements of the total thermal energy in clusters, Y, strongly 
depend on the projection axis through the cluster. Thus, mis- 
alignment between the semi-analytic baryon distortion and the 
"true" distribution may cause biases when using semi-analytic 
models to, e.g., tie weak-lensing and SZ observations together 
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Figure 10. The dependence on mass and the redshift evolution of the cumulative ellipticity profile as a function of r/RoQQ. Left: Shown is the ellipticity profile at 
z = for various mass bins. Bottom left: Shown are the relative differences in ellipticity to the lowest mass bin (1.1 X IO'^^Mq < A/200 < 1-7 X IO'^'Mq ). Over 
this mass range, the cluster ellipticities show a noticeable but not substantial mass dependence within S500. in contrast to the stronger dependence on Pkin/^th- Right: 
Shown is the ellipticity profile for various redshift bins. The horizontal purple and pink error bars have the same meaning as in Fig. [T] Bottom light: Shown is the 
relative difference of ellipticity at a given redshift to z = 0. The redshift evolution of the ellipticity (especially at large radii, r > Sjoo) is driven by the larger amount 
of substructures at higher redshifts due to the increased mass accretion rate of group/cluster halos at these redshifts. The pressure-weighted ellipticities track the 
density-weighted ellipticities well and show the same trends with redshift. 
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Figure 11. The weighted median angles between the DM major axis and gas 
(red) and pressure (green) axes at each radius for all simulated physics models: 
AGN feedback (solid), radiative cooling (short-dashed), and shock heating- 
only (long-dashed). The 25* and 75"' percentile values are shown for the gas 
density in the AGN feedback model (dotted). On average the gas and pressure 
axes are misaligned by 20 to 30 degrees to the DM principle axis, independent 
of our simulated physics models. However, both simulations with radiative 
cooling show more misalignment in the inner regions than the non-radiative 
simulations. The hght colors and lines represent the region where the average 
cluster shape is close to spherical (f/a|DM > 0.75) such that the major axes are 
not well defined and their angles are approaching a random distribution. Note 
that we have weighted the average angles by l-c/a to down- weight the angles 
from the spherical ICM shapes and the cluster interiors. The horizontal purple 
and pink error bars have the same meaning as in Fig. [T] 



5. SZ SCALING RELATION 

In this section we explore the impact of AGN feedback, clus- 
ter shapes and kinetic pressure support on the SZ-flux-to-mass 
relation, Y-M, using our large sample of clusters. We compute 
the SZ flux for all clusters for both, spherical boundaries and 
cylindrical apertures (Fsph and Fcyi)- For the cylindrical aper- 
ture calculations the total fluxes are computed along each axis 
of the moment-of-inertia frame, measured at /?200, and addi- 
tionally along each axis of another randomly-oriented frame. 
We choose the line of sight boundaries for the cylindrical in- 
tegrations to be three times the radius of the aperture. This 
procedure enables quantifying the importance of substructure, 
which we have already shown in Sections[3]and|4]to be signifi- 
cant at radii beyond 7? 200- From the calculated Fa values we fit 
an average scaling relation. 



Fa = 10^ 



Ma 



3 X W^h-lMQ 



/^vdMpc^ 



(10) 



where A and B are the fit parameters for the slope and normal- 
ization, respectively. We weight each cluster by its Fa when 
fitting for A and B to keep the low -mass clusters from com- 
pletely dominating the fit. 



5.1. Self-similar Y-M scaling relation 

We review the expectations for F in the idealized case of 
a cluster in virial equilibrium to help understand how possi- 
ble deviations from the self-similar Y-M relation and the scat- 
ter about it may arise. Starting with Eq. (O, which has been 
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Figure 12. The Y-M scaling relation for the AGN feedb ack simulations 
compared to recent X-ray results from 'Amaud et al.l iSOTOl) and SZ results 
from ACT (Marriage et al. 2010), SPT ( Andersson et al. 2010), and Planck 
jPlanck C ollaboration et al. 2011a). We have applied the 15% correction to 
the X-ray A/hse from lKravtsov et al.l (200611 . 



rewritten as, 



Y = 



(7-1)ctt 

dV Pe = ^ JCe M -Egas , 



(11) 



where is the electron fraction defined as the ratio of electron 
and hydrogen number densities = n^/riii = (Xh+ 1)/(2Xh) = 
1 . 158, 7 = 5/3 is the adiabatic index, ^ = 4/(3Xh+ 1 +AX^x^) = 
0.588 denotes the mean molecular weight for a fully ionized 
medium of primordial abundance, and we assume equilibrium 
between the electron and ion temperatures. Next, we define 
the ch aracteristic temperature of the halo fKomatsu & Sehald 
l2002l) as 



^^200 ■ 



3^200 3 



[10G//oM2oo£(z)] 



2/3 



(12) 



so we can write the total thermal energy of the halo with 
Eq. (fT2l i as 



Egas - 2 -^ga-s ^^200 - ( 1 



-J*)jbjc 



2R 



200 



= (i-/*)./b./c ^ [mop,Az)]'^'Mll^y 



(13) 



Here /» < M^,/M\, is the stellar mass fraction within the halo 
and fc is the correction factor for the fraction of missing 
baryons at a given overdensity. Then we insert Eq. ( fT3] ) into 
Eq. (fTTT i to get the integrated Compton-y parameter within 

^200, 



XeXiin(l-f*)fbfcG 



1/3 



M: 



5/3 



200 



M200 



5/3 



f7b 0.25 
0.043 



(14) 



For Eq. ( fT4l i. we set = 0, /c = 0.93 (as calculated from our 
shock heating simulations at /?20o) and adopted the cosmolog- 
ical parameters of our simulation. This simple analytical ex- 
pression for the Y-M scaling relation allows one to explore the 
assumptions underlying its derivation. More specifically, we 
test the assumptions of spherical gravitational potential, zero 



non-thermal pressure support, and constant fi, (and for simu- 
lation with star formation, constant /*) at Ra, independent of 
cluster mass. 

5.2. Comparison to data 

In Figure [121 we compare Y^oh for our simulated clusters to 
the X-ray results fromjArnaud et al. (2010), and the SZ re- 
sults from ACT ("Marriage et al. 2010), SPT dAndersson et alj 
120101) . and Planck (Planck Collaborati on et al.1 1201 lal) . ~We 
adopt the 15% correction to the X-ray Mhse estimates from 
Kravtsov et al. (2006) which is valid for the respective obser- 
vational sample selection criterion. Our Fsph-M relation with 
AGN feedback is consistent with the current data from X-ray 
and SZ observations. However, at group scales, our simulations 
slightly overpredict the SZ flux due to the too high gas frac- 
tions, /gas, in our simulations compared to X-ray observations 
(cf. BBPS4). Potentially our simulations are missing some of 
the re levant physics that governs /g^s (see, e.g.. lPfrommer et al] 
1201 Ih or underestimate the action of AGN feedback on these 
mass scales. 

The Fsph reported by SZ surveys for known clusters use an X- 
ray-derived estimate of the aperture size. This is useful because 
the cluster radii are typically poorly measured in SZ, and so the 
X-ray aperture fixes the SZ measurement along the otherwise 
degenerate aperture flux/aperture radius relation. However, this 
prior introduces correlations between the X-ray and SZ obser- 
vations, which makes comparisons between these observations 
difficult to interpret. 

5.3. Physics dependence of the Y-M relation 

In Figure[T3]we show the dependence of the Y-M relation on 
our three simulated physics models, i.e., shock heating, radia- 
tive cooling and star formation, and AGN feedback. The stark 
differences between the shock heating simulation and the two 
radiative simulation models arise from the loss of baryons in 
the ICM to star formation. The radiative cooling simulations 
show a constant normalization offset of ^ 20%, which nearly 
matches the values for these simulations. In Table|2]we show 
that the self-similar expectation of Eq. ( fT4b almost completely 
captures the cluster thermodynamics in our simulations when 
integrated over cluster-sized apertures. Including more physi- 
cally motivated sub-grid models in the simulations, we find that 
both, the shock heating and radiative cooling slopes are consis- 
tent with this self-similar derivation for the Y-M relation, while 
the AGN feedback simulations have a steeper, mass-dependent 
slope. This break from self-similarity in the AGN simulations 
arises from the suppression of star formation in the higher mass 
clusters and a feedback-induced deficit of gas inside the lower 
mass clusters. Over the redshift ranges we explore (z = to 
z = 1.5) and for all simulated physics models, the Y-M scaling 
relation normalization changes as predicted by self-similar evo- 
lution and the slopes remain essentially constant (cf. Fig.[T4b. 
So, the Y-M relations from AGN simulations are different at 
z = 0, but evolve as predicted by self-similar evolution. This 
result is independent of the two aperture sizes chosen, which 
correspond to over-densities of 200 and 500 times the critical 
density (Fig. [T4l l. As we have repeatedly found in previous 
sections, the clusters interior to the radii Rsm and even 7?2oo 
are relatively well behaved, with only modest impact of cluster 
ellipticities and kinetic pressure. 

To quantify the scatter, we compute the relative devia- 
tion of each cluster from the mean relation, SY^ /Y = (1a - 
J'A,fit)/5'A,fit, and then fit this distribution with a Gaussian prob- 
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Table 2 

Y-M scaling relation fits for different simulated physics, sub-sampling in kinetic-to-thermal energy and ellipticity (of the density and pressure distribution), and 

along different projected axes yielding Fcyl- 



A = 200 


1 = 


z = 0.5 


z= 1 




B A cry 


B A (Jy 


B A cry 










Theory, Eq. {14) 
Shock heating 

R nrii nti vf rrinlin o 

AGN feedback 


-4.88 1.67 
-4.87 ±0.01 1.64 ± 0.03 0.115 ±0.001 
-4 Q4 -1- 01 1 67 + 03 1 1 9 + 001 
-4.92 ± 0.01 1.71 ±0.03 0.133 ±0.001 


-4.81 1.67 
-4.81 ±0.02 1.63 ±0.04 0.116 ±0.001 
.4 00 + n no 1 67 + 05 121+ 002 
-4.87 + 0.02 1.72 ±0.05 0.136 ±0.002 


-4.74 1.67 
-4.76 + 0.05 1.61 +0.07 0.111 +0.002 
.4 oT _i_ n 05 1 66 + 09 121+0 002 
-4.82 + 0.05 1.73 + 0.09 0.147 + 0.003 


K/U" 








Lower 3"* 
MiHHIp 3'"^ 

iVllUUlC J 

Upper 3'^'^ 


-4.88 ±0.04 1.73 ± 0.07 0.115 ±0.002 
.4 Q2 + n 02 1 73 + 05 116 + 002 
-4.94 ±0.02 1.72 ± 0.04 0.123 ±0.002 


-4.83 ± 0.06 1.73 ±0.11 0.122 ±0.003 
.4 07 + n 04 1 7:! + n OR 121+0 003 

-4.89 + 0.03 1.73 + 0.07 0.122 + 0.003 


-4.78 + 0. 14 1 .76 + 0.2 1 0. 1 37 + 0.004 
.4 09 + 07 1 73 + 1 3 138 + 005 
-4.85 + 0.09 1.72 + 0.15 0.141 +0.005 


c/a (gas)" 








Lower 3"* 
Middle 3'^'' 
Upper 3"* 


-4.94 ±0.02 1.72 ±0.05 0.131 ±0.002 
.4 go -1- 02 1 70 + 05 1 29 + 002 
-4.90 ±0.02 1.72 ± 0.05 0.119 ±0.002 


-4.89 + 0.04 1.72 + 0.08 0.133 + 0.003 
.4 o/s + 05 1 73 + 09 133 + 003 
-4.86 + 0.04 1.72 + 0.08 0.123 + 0.003 


-4.84 + 0.09 1.73 + 0.15 0.147 + 0.005 
.402 + 010 1 72 + 016 01 36 + 004 
-4.80 + 0.09 1.73 + 0.15 0.138 + 0.004 


c/a (pressure)" 








Lower 3''' 
Middle 3'''' 
Upper 3"^ 


-4.94 ±0.02 1.71 ± 0.05 0.132 ±0.002 
-4.92 ±0.03 1.71 ±0.05 0.129 ±0.002 
-4.90 ±0.03 1.72 ±0.05 0.122 ±0.002 


-4.88 + 0.04 1 .73 + 0.08 0. 1 34 + 0.003 
-4.87 + 0.04 1.72 + 0.08 0.128 + 0.003 
-4.86 + 0.05 1 .7 1 + 0.09 0. 129 + 0.003 


-4.84 + 0.08 1.73 + 0.14 0.152 + 0.005 
-4.82 + 0.08 1.72 + 0.14 0.134 + 0.004 
-4.78 + 0.12 1.76 + 0.19 0.137 + 0.004 


b/a (gas)° 








Lower 3''' 
Middle 3''' 
Upper 3"^ 


-4.93 ±0.02 1.70 ±0.04 0.138 ±0.007 
-4.92 ±0.02 1.71 ± 0.05 0.130 ±0.006 
-4.91 ±0.02 1.72 ±0.05 0.128 ±0.002 


-4.88 + 0.04 1.72 + 0.07 0.140 + 0.010 
-4.87 + 0.04 1 .72 + 0.08 0. 1 37 + 0.009 
-4.86 + 0.04 1 .72 + 0.07 0. 1 32 + 0.003 


-4.83 + 0.07 1.73 + 0.12 0.149 + 0.016 
-4.83 + 0.11 1.72 + 0.18 0.149 + 0.014 
-4.81 +0.09 1.73 + 0.14 0.136 + 0.004 


Ycy\ rotated" 








Minor axis 
Middle axis 
Major axis 


-4.87 ±0.01 1.69 ±0.03 0.126 ±0.001 
-4.87 ±0.01 1.69 ± 0.03 0.128 ±0.001 
-4.84 ±0.01 1.68 ±0.03 0.137 ±0.002 


-4.83 + 0.02 1 .69 + 0.05 0. 1 36 + 0.002 
-4.82 + 0.02 1 .69 + 0.05 0. 1 33 + 0.002 
-4.79 + 0.02 1.67 + 0.05 0.152 + 0.002 


-4.79 + 0.05 1.70 + 0.09 0.144 + 0.003 
-4.78 + 0.05 1.69 + 0.09 0.147 + 0.003 
-4.74 + 0.05 1.68 + 0.09 0.167 + 0.003 


Fcyi random" 








axis 1 
axis 2 
axis 3 


-4.86 ±0.01 1.69 ±0.03 0.129 ±0.001 
-4.86 ±0.01 1.69 ±0.03 0.129 ±0.001 
-4.86 ±0.01 1.69 ± 0.03 0.130 ±0.001 


-4.82 + 0.02 1.69 + 0.05 0.140 + 0.002 
-4.82 + 0.02 1 .69 + 0.05 0. 1 38 + 0.002 
-4.82 + 0.02 1.69 + 0.05 0.141 +0.002 


-4.77 + 0.05 1.70 + 0.09 0.146 + 0.003 
-4.78 + 0.05 1 .68 + 0.09 0. 150 + 0.003 
-4.77 + 0.05 1 .70 + 0.09 0. 152 + 0.003 



" For fits to all sub-samples/projections, we use our AGN feedback model. Fit parameters are defined in Eq. HOi . 
* For the h/a sub-sampling of Y^yi, we chose random axis 1. 




Figure 13. The normalization, slope and scatter of the Y-M scaling relations all depend on the simulated physics. Left: The Y-M scaling relations at ^ = for all 
simulated physics models: shock heating (green), radiative cooling (blue), and AGN feedback (red). The )'-axis has been scaled by M^^' to highhght the deviations 
from self-similarity. Right: The probability distributions for the relative deviation, 5Yj\/Y, with respect to the best fits for all three physics models. We also show 
Gaussian fits (dotted lines) and include Poisson deviations for the AGN feedback simulations (grey band). We find that the AGN feedback simulations have the 
largest scatter and a steeper slope compared to the other simulations. 
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Figure 14. All simulated Y-M scaling relations evolve self-similarly with redshift according to Eq. )14t . We show the Y-M scaling relation fits for the normalization, 
B, (left panel) and slope. A, (right panel) as a function of redshift and for two different cluster masses M200 and -M500, and compare those to the self similar prediction 
for M200 (dotted black). The Y-M relation of AGN feedback simulations has a different slope, but shows no anomalous redshift evolution relative to self-similar 
evolution. 




Figure 15. Rotating the clusters into their major, middle and minor axes and calculating projected (cylindrical) Y-M relations shows the effect of infalling substruc- 
ture. Left: The cylindrical Y-M scaling relations from the AGN simulations for clusters that have been rotated into their major, middle, and minor axes defined by 
computing the (3D) moment-of-inertia tensor within Rim- Right: The probability distributions for the scatter, SY/^, relative to the best fits for all three distributions, 
each representing a distinctive rotation as well as the spherical distribution (black). We include the Gaussian fits (dotted lines) and the Poisson deviations for the 
major axis rotation (grey band). Rotating the clusters such that integration happens along the major axis increases the total Y values, while further distorting and 
increasing the scatter (due to the large cluster-to-cluster variance in the infall regions). The Y^yi values are integrated along the given axis from — 3R200 to 3/?20o; 
hence for any given cluster F^yi > Fsph ■ 
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Figure 16. Sub-sampling the Y-M relations by the kinetic-to-thermal energy ratio (K/U) for the AGN simulations. Left: The Y-M scaling relation for the three 
K/U sub-samples, upper 3"^ (red), middle 3"^ (green) and lower 3'^'' (blue), with the corresponding slope fitted to those points. The v-axis has been scaled by M^^^ 
to highlight the deviations from self-similarity. Right: The probability distributions for the relative deviation, 5Y/\/Y, with respect to the best fits for the three 
sub-samples and the total distribution (black), including the Gaussian fits (dotted lines) and the Poisson deviations for the upper 3"* sub-sample (grey band). The 
sub-sample of K/U with the largest kinetic pressure support (upper 3"*) shows systematically lower total Y values for a given mass as well as larger scatter, while the 
lower K/U sub-sample has the lowest scatter of ~ 11%. This is expected because a larger kinetic pressure contribution implies a lower thermal pressure and hence 
decreases Y. 



ability distribution function (PDF), 
G(SYA/Y)=Aoexp 



-(SYa/Y) 



21 



(15) 



Here the parameter Ao is the normaHzation and uy is the vari- 
ance, which we will refer to as the scatter Here we have chosen 
to model the va riation about t he me an as a Gaussian, while pre- 
vious work by IStanek et alj (120101) showed that a log-normal 
distribution is also a reasonable description of the scatter In 
Appendix iDl we show that within the (Poisson) uncertainties, 
the scatter is clearly Gaussian distributed and only approxi- 
mately log-normal. Forcing a log-normal distribution intro- 
duces higher-order moments such as skewness and kurtosis as 
can be seen by the tails in the distributions and their asymmet- 
ric shapes. 

We find that the scatter, cry, for the entire sample of clusters 
is between 11 % and 13 % (cf. Fig. [TSjand Table |2li, which is 
consistent with previous work {Nag ai 2006; Shaw et al. 2008; 
IStanek et al"]l2010l: lYanget all 12010) . In the simplest simula- 
tions with only shock heating the source for this scatter in the 
Y-M relation has been proposed to arise from the formation 
time, the co ncentration, and the dynamical state of the cluster 
(lYangetal.l f2010:). As our simulations include more sub-grid 
physics models the scatter increase from ~ 11% to ^13% at 
z = and changes further from ~ ll%to~ 15%atz = 1. Of the 
three different physics models, the simulations with AGN feed- 
back gives the l argest scatte r, which is consistent with semi- 
analytic results (Shaw et al. 2008). This model for AGN feed- 
back is self -regulated (Battaglia et al. 20 10) and injects ^2/3 
of the energy before z = 1 when the average cluster mass is sig- 
nificantly smaller and the associated potentials are shallower so 
that a fixed energy injection by AGNs may in principle have a 
stronger impacQ Thus, the increased scatter in the Y-M rela- 
tion from the AGN feedback simulations compared to the sim- 

^ Similar results were found bv lMcCarthv et al] 1201 II) in simulations with 
a more detailed feedback prescription. 



ulations without feedback is a result of the energy injection, 
which heats and disturbs the ICM. This statement is in accor- 
dance with previous results from Battagl ia et al.l (|2010), where 
they showed the impact of AGN feedback on the pressure pro- 
files of clusters and found that simulations with feedback had a 
shallower asymptotic pressure profile slopes than those without 
feedback. Thus, the intermittent nature of energy injection into 
the group system early-on results in a larger scatter in the Y-M 
relation compared to simulations without energetic feedback. 

5.4. Cylindrical Apertures 

For pointed SZ observations of clusters and SZ surveys, a 
natura l, model-independent observable is the proj ected flux, 
Feyi (iMroczkowski et all l2009t ISavers et all l20Tll) . We find 
Jcyi > J'sph in all cases, whether we chose the projection along 
a principal or a random axis. This is due to the assumed ex- 
tension along the line-of-sight integration which we choose to 
be three times the aperture radius; in observations, structure 
beyond this scale may additionally contribute in some cases. 
In fact, a projection integral out to 37?2oo decreases the Y-M 
slope for the AGN feedback simulations such that it becomes 
consistent with the self-similar slope (cf. Table |2]l. We find no 
difference between the random 2D projections and the integra- 
tion along the middle or minor axes with respect to the nor- 
malization and slope (cf. Fig. [15] and Table |2]i. The scatter for 
the random 2D projections is marginally larger than the projec- 
tions along middle and minor axes. Our results show that the 
integration along the major axis yields dramatically different 
results, both, for the normalization and scatter in comparison 
to projections along the other axes. This has its origin in the 
more extended tails of the PDF (cf. Fig [TsT l. The normaliza- 
tion and scatter between the major axis and the other axes in- 
crease by ^ 7% and ^ 9%, respectively. At a higher redshifts, 
these differences are amplified and we find a ^ 12% increase 
in the normalization and an increase in scatter by ^ 16%. This 
indicates that substructure is preferentially aligned with major 
axis and that substructure heavily influences the result from the 
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Figure 17. Sub-sampling the Y-M relations by the gas c/a axis ratio (upper panels) gas h/a axis ratio (lower panels) for the AGN simulations. Left: The Y-M 
scaling relation for the three c/ci sub-samples, upper 3"* (red), middle 3'^'' (green) and lower 3'^'' (blue), with the corresponding slope fitted to those points. The y-axis 
has been scaled by M^/^ to highlight the deviations from self- similarity. Right: The probability distributions for the relative deviation, (57^/7, with respect to the 
best fits for the three sub-samples and the total distribution (black), including the Gaussian fits (dotted lines) and the Poisson deviations for the upper 3"* sub-sample 
(grey band). The sub-sample of c/a containing the lowest values (largest ellipticities) shows systematically lower total Y values for a given mass and larger scatter, 
while the more spherical high c/a sub-sample shows a lower scatter of ~ 12%. Additionally we find that the pressure c/a axis ratio sub-sample has similar results 
(cf.Table|2). The sub-sample oib/a containing the lowest values (largest projected ellipticities) shows larger scatter than the sub-sample with the highest values by 
~ 23%. There are significant changes in the Y-M relation at z = 0, however, smaller in comparison to the c/a sub-sampling. 



moment-of-inertia tensor beyond /?200- 

5.5. Toward a fundamental plane of Y-M 

After quantifying the scatter of the entire sample, we aim at 
understanding its origin. This may enable us to either construct 
a linear combination of physically motivated observables that 
minimizes the scatter or to employ sub-sampling of the full dis- 
tribution according to some parameter so that the resulting dis- 
tribution exhibits a smaller intrinsic scatter an d potentially al - 
lows for tighter cosmological constraints (e.g., Afshordi 2008). 

In the previous sections we explored the average radial trends 
for kinetic pressure support from bulk motions and gas den- 
sity/pressure shapes of the ICM. Utilizing this information, we 
rank order clusters according to their kinetic pressure support 
and intrinsic shape information. We follow the same fitting 
procedure as above for subsets of the lower 3'^'', middle 3"^, and 
upper 3"^ of the correspondingly sorted distributions in order 



to demonstrate the impact of kinetic pressure support and as- 
phericity on the Y-M relation fits and scatter. For the rest of this 
section we concentrate our analysis on the Y-M relations of the 
AGN feedback simulations, since they show the largest scatter 
(this will provide an upper limit on the scatter) and are most 
likely our best representation of "real" clusters in comparison 
to the other simulated physics models. We compute the ratio of 
kinetic-to-thermal energy, K/U , within radial bins and use this 
ratio as a measure of dynamical state for the galaxy clusters. 
We define the internal kinetic energy, K, and thermal energy, 
f/, of a cluster as 



g«r)^^ ^"^7^'""'' , (16) 
U{<r)^Y.^-^^^f^^ (17) 
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Figure 18. The correlations between the K/U ratio and 1 —c/a. Here the red 
points represent each cluster in the simulations and the blue crosses are average 
quantities. The linear correlation coefficient is 0.575 

where m and p are the gas mass and the SPH density, respec- 
tively for all particles / less than radius r. The ratio K/U is the 
volume integrated analog of the ratio i\in/Ah shown in §[3] and 
hence is also an indicator of formation history and substruc- 
ture. For the sub-sample with the highest ratio of K/U, we find 
a smaller normalization (cf. Fig. [16] and Table |2]l where the dif- 
ference between this upper 3"^ of the distribution and lower 3"* 
is ~ 15%. Here some of the thermal pressure support has been 
compensated for by kinetic pressure support resulting in lower 
integrated thermal electron pressure, thus, lowering F-values. 
More massive clusters are typically in the high K/U sample 
rather then the other two samples. We find that the sub-sample 
with the smallest K/U values shows the lowest scatter, ~ 1 1% 
for the AGN feedback simulations. Further sub-sampling of the 
smallest K/U values (e.g., the lowest 6*) does not decrease the 
scatter, which is limited to 11%. 

We also sort our cluster sample by the ratio of minor to ma- 
jor axis for both 3D and projected 2D, c/a and b/c^ respec- 
tively, as defined in Sec. H) Following the same procedure as 
for the K/U sub-sample and restricting ourselves to the AGN 
feedback simulations, we find that splitting the clusters up by 
ellipticity, c/a, gives similar results in comparison to K/U- 
splitting. The galaxy clusters with smaller ellipticities have 
larger total Y values and less scatter, while the more triaxial 
clusters have lower total Y and large scatter (cf. Fig. [TT] ). These 
trends are reflected in the fit parameters of the sub-sample Y- 
M relation shown in Table |2l where the differences between 
the upper 3"* and lower 3'^'' sub-samples normalization param- 
eters is ^ 12%. Additionally, we found that using the pressure 
shapes instead of the gas shapes yield almost identical results 
(cf. Table |2]i. We find that sub-sampling the clusters with the 
b/a statistics has the greatest impact on the scatter The sub- 
sample of clusters that appear to be elongated in the plane of 
the sky have larger scatter than the more spherical clusters (cf. 
Fig. [l7l and Table |2]i. Also, the b/a sub-sampling causes a 
similar bias in the Y-M relation in comparison to the c/a sub- 
sampling, but not as significant. 

The results from sub-sampling clusters according K/U and 
c/a indicate that there are correlations between these physical 

^ We show one projection for the b/a sub-sampling. The other projections 
yield similar results. 



properties and the scatter in the Y-M relation. In Figure [18] 
we show that larger K/U ratios correlate with larger 1- c/a, 
i.e. larger triaxiality, with a linear coefficient value, = 0.58. 
These correlations between kinetic pressure support and ellip- 
ticity are the result of the growth of structure being hierarchical. 
This supports the argument that kinetic pressure support, ellip- 
ticity and sub-structure are all tracers of the dynamical state 
and the formation history of galaxy clusters, which is the ulti- 
mate cause of the intrinsi c scatter of the Y- M re lation. Simi- 
lar res ults were found by Ra sia et aTl (1201 Ih and iKrause et alj 
(1201 ll) Previous work by Yan g et alJ (1201 Oh found mass trends 
in the measured scatter, which is consistent with our findings 
after extrapolating their lower mass range to our larger masses. 
However, their conclusion is different from ours, since they 
claim that the scatter is most sensitive to the DM concentration; 
a finding that may partially be due to the insufficient resolution 
in their simulations. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we demonstrate that the spatial distribution of 
the ICM, kinetic pressure support from bulk motions, and self- 
regulated thermal energy feedback in clusters cores (that we 
refer to as AGN feedback) all play very important roles for 
the thermal properties of galaxy clusters. In particular, the ob- 
servables for large SZ cluster surveys, such as ACT, SPT and 
Planck, will be modified by these processes. Below we high- 
light and expand on our main results. 

Non-thermal pressure support and cluster shapes: The con- 
tribution to the overall pressure support in clusters from bulk 
motions, Pyn, increases substantially for larger radii and is a 
strong function of both, cluster mass and redshift. Includ- 
ing AGN feedback marginally decreases Pkin/Ah in compari- 
son to the other (more simplified) simulation models, namely 
our shock heating-only model and that which additionally in- 
cludes radiative cooling, star formation, supernova feedback, 
and CRs. However, the difference is not substantial enough to 
be statistically inconsistent with the variance around the me- 
dian of Pkin/Ah- The mass dependence and redshift evolution 
of fldn/Ah is governed by P^n and a direct result of the hi- 
erarchical growth of structure. Semi-analytic approaches are 
just beginning to model P^n. The full dependence on radius, 
mass and redshift of this component is, by definition, self- 
consistently included in hydrodynamic simulations. 

We find that the distribution of gas density and pressure are 
weak functions of the simulated physics models within P200 
(excluding the cluster core) and that AGN feedback mildly 
modifies the average gas shapes. The cluster mass dependence 
of the ellipticity is more moderate in comparison to Pkin/Ah- 
The ellipticity is small within P500 with little redshift evolution. 
In combination with the comparably small non-thermal pres- 
sure support at these scales (which rises dramatically beyond 
this characteristic radius), the small clumping factor measured 
in our simulations (cf. BBPS4), and the small modification 
of our simulated cluster physics at these radii (in particular of 
our implementation of AGN feedback), this result is reassur- 
ing for X-ray observations of clusters which use ^500 to char- 
acterize clusters with high-quality Chandra and XMM Newton 
observations. Hence, our analysis theoretically supports this 
choice of radius (which was initially motivated by the simu- 
lations in Evrard et al. ( 1996)) and justifies some of the main 
assumptions such as spherical symmetry and an almost radius- 
independent hydrostatic mass bias of ^ 20-25% when using a 
fair sample of clusters without morphological selection which 
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may be applicable for the future eROSITA sample. 

We find substantial redshift evolution in different dynamical 
quantities, e.g., Pkin/Ah, the velocity anisotropy, and anisotropy 
parameters such as ellipticities. This is in particular the case for 
the changes in the power-law behaviors of the radial profile of 
these quantities such as the sudden break in ellipticities which 
moves to smaller radii as the redshift increases (when scaled 
to 7?2oo)- The break and the more pronounced ellipticities and 
7\in/Ah outside a characteristic radius are a direct result of in- 
creased level of substructure predicted by hierarchical struc- 
ture formation and the associated higher mass accretion rate at 
higher redshift. We explicitly show (in the Appendix |C]l that 
most of this redshift evolution is somewhat artificial and can 
be absorbed in a re-definition of the virial radius: scaling with 
the radius that contains a mean density of 200 times the aver- 
age mass density rather than the critical density of the universe 
considerably weakens the observed trends with redshift. This 
also suggests a physical definition of the virial radius in terms 
of dynamical quantities (that, however, remain poorly defined 
observationally), e.g., the equipartition radius of thermal and 
kinetic pressure, the region where the velocity anisotropy be- 
comes strongly radial, or the radius at which the ellipticity or 
substructure level increases dramatically. These seemingly dif- 
ferent criteria all select a rather similar radius around ^200,m; 
almost independent of redshift. 

On scales > /?5oo, stacking analyses of projected SZ cluster 
images can be done with data from SZ experiments such as 
ACT and SPT, and, provided there is a suitable sample size, 
one may be able to detect projected gas pressure shapes, po- 
tentially even in bins of redshift. The results on the randomly 
projected 2D axis ratios represent the theoretical expectations. 
Any statistics from the intrinsic 3D distribution is highly corre- 
lated with the projected 2D distribution; we find that the (more 
elliptical) intrinsic cluster shapes can on average be inferred 
from their projected analogues by applying a ^ 5 - 10% cor- 
rection on the ellipticity. Another interesting outcome from 
our shape analysis is that there is no direct and simple map- 
ping of shapes and alignments for DM spatial distribution to 
the gas and pressure distributions possible mostly due to the 
difference in substructure distribution and dissipational nature 
of the gas. This result is troublesome for semi-analytic mod- 
els which use dissipationless simulations as a template to solve 
for the gas distributions and pressure shapes. Such a method 
will produce additional triaxiality and misalignment for such 
a semi-analytical model of the ICM. The overall magnitude of 
the shape is reconciled by using the gravitational potenti al (e.g., 
lOstriker et al]|200l iBode et al.ll2009l: iTrac et alJl20Tl which 
has been shown to be less triaxial (lLauetalJl20'TW " than the 
DM. However, providing an algorithm to re-alignment these 
pseudo gas distributions is a non-trivial task. 

Y-M scaling relations: Our simulations are in good agree- 
ment with the current Y-M scaling relations from both X-ray 
observations and SZ surveys. However, to properly predict the 
Y-M scaling relations for an SZ experiment such as ACT, SPT 
or Planck without any prior knowledge of cluster masses, care- 
ful mock observations are needed. Those would have to in- 
clude a simulation of the CMB sky with associated experiment 
noise and adopt the relevant cluster selection pipelines for the 
given experiment that employs the same cluster profile used for 
matched filtering in order to include all the systematics and po- 
tential biases that are intrinsic to the data analysis, e.g.. X-ray 
priors on the aperture size. 

We find that the inclusion of AGN feedback causes a de- 
viation from the predictions of self-similar evolution for both 



the normalization and slope of the Y-M relation (as measured 
within 7?2()o)- However, we recover the self-similar slope again 
in our projected Y-M scaling relations (where we integrate 
along a cylinder of half-height 3/?2()o), suggesting that AGN 
feedback pushes a fraction of its gas beyond the virial radius 
and a larger aperture/projection radius is able to recover the 
thermal energy from this larger reservoir of gas. 

Including AGN feedback also increases scatter in the Y-M 
relation compared to simulations that include shock heating 
alone, from ^11% to ^13%. Interestingly, sorting the clus- 
ters into sub-samples of K/U and c/a will reduce this scat- 
ter; e.g., K/U sub-sampling reduces the scatter from ~ 12% 
to ^ 11%. We find that sub-sampling introduces only a small 
(predictable) bias in the normalizations on the order of a few 
percent. This suggests that observational proxies for the dy- 
namical state and ellipticities may be used to construct a funda- 
mental plane of the Y-M relation. The scatter ultimately orig- 
inates from the merging history with its redshift and mass de- 
pendent accretion rates; those determine the non-thermal pres- 
sure support, the level of substructure, and the ellipticity. While 
sub-sampling on one of these secondary tracers may decrease 
the scatter, it is unlikely to decrease much more if more tracers 
are used (as they probe the same underlying process, albeit with 
a different weighting). Conversely, our band-splitting analysis 
on the Y-M relations suggests that large outliers from the mean 
relations would be interesting candidates for follow up with 
high resolution SZ observations, since they are more likely to 
have larger kinetic pressure support and ellipticities. 

A fundamental point to take away is that all results at larger 
radii (> /?2()()) for the kinetic pressure support and ICM shapes 
are dominated by substructure. We also see the impact of sub- 
structure on the cylindrical Y-M scaling relation when integrat- 
ing along the major axis with which substructure is preferen- 
tially aligned. Quantifying substructure statistically is difficult 
because of the problem of double-counting: the large volume 
contained within the radius that contains 95% of the total SZ 
flux, 4/?200, necessarily leads to overlapping volumes of neigh- 
boring clusters, especially at high-redshift. Thus, this property 
remains challenging to model phenomenologically or analyti- 
cally. 

As discussed previ ously in the literature (e.g. lBattagha et all 
I20T0I: ISunetalJl20f l) SZ galaxy cluster may provide further 
insight into the interesting astrophysics associated with the 
ICM of clusters. This however may significantly complicate 
cosmological analyses in producing competitive constraints. 
However, these are exciting prospects for studies of feedback 
and other energy injection processes within clusters especially 
at higher redshift since the selection function of SZ cluster sur- 
veys probes clusters which populate the massive and high red- 
shift end of the distribution. 
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Figure 19. The kinetic pressure-to-total pressure is weakly mass-dependent, 
Pjjin/Ptot oc ^200 indicated by the scaUng of the y-axis. Shown is the median 
of ^kiii/^tot as a function of radius for the AGN feedback simulations for vari- 
ous mass bins with the 25* and 75"' percentile values illustrated by the dotted 
lines for the lowest mass bin at z = 0. For comparison, we also show the model 
for /\in/^tot by Shaw et al. (2010), which has been fit to match AMR simula- 
tions (dash-dotted). We illustrate the 1 and 2 a contributions to F/\ centered 
on the median for the feedback simulation by horizontal purple and pink error 
bars. Therefore, ignoring this mass dependence results in a 60 % difference in 
this ratio for an order of magnitude change in the cluster mass. The median of 

PvinjPta scales as 'Wi™, which results in a larger difference. 




Figure 20. Shown is the relative difference between axis ratios with and with- 
out the r"^ weighting for the gas-density (red line), DM-density (blue line) and 
gas-pressure (green line) weightings. Additionally including the r"'^ weighting 
in the definition of the moment-of-inertia tensor down-weights the contribution 
at larger radii by < 15%, thus reducing the effect of substructure. 

APPENDIX 
A. FITTING FUNCTION FOR Agn/^Vot 

In Section [3] we show that the ratio /"kin/Ah is a function 
of mass. However, the prev ious empirical fitting function for 
fldn/Aot (IShaw et al.l 120101) does not include a mass depen- 
dence. 
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where a{z) = ao(l +z)^ for low redshifts (z < 1) and the fit 
parameters are ao = 0. 18 ± 0.06, /3 = 0.5, «nt = 0.8 ± 0.25, and 
by construction, »m = 0. In Fig. [19] we compare the fitting 
function for Eq. (lAll) and i\in//^ot, split by different mass bin 

which have been scaled by Mjqq, i.e. hm = 1/5, that minimizes 
our x^. We chose a normalization of 3 x 1O''*M0 to match 
the fitting function of IShaw et al.l (120101) . Thus, the median 
difference between Pkin/Aotof a 1O'^M0 and a lO'^^Mo cluster 
is ^ 60%. The ratio in Eq. (lAll l is similar to what is shown Fig. 
[1] however, Pkin/Ah depends more sensitively on mass. We find 

that the mass dependence for this ratio amounts to M- 
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DOWN-WEIGHTING THE SUBSTRUCTURE IN THE 
MOMENT-OF-INERTIA TENSOR 



(Al) 



For both the gas density- and pressure-weighting of the 
moment-of-inertia tensor, the inclusion of an additional x~^- 
weighting has a relatively minor influence on cluster shapes 
(cf. Fig. l20b and we do not see large differences in the axis 
ratios at the larger radii. The x~^-weighting does lessen the in- 
fluence of substructure which we have seen to be important at 
radii beyond /?200i but it does not remove it or isolate its signal. 
This would be a non-trivial ta sk for any stack ing analysis as it 
was recently suggested by Zemp et al.l (1201 lb . 

C. CLUSTERS IN VELOCITY SPACE AND A DYNAMICAL 
RADIUS DEFINITION 

The radial trends over redshift seen in Figures [Tl and [TOl call 
for re-examination of the choice for the working definition of 
radius, which is directly r elated to the definition of the clus- 
ter mass (see IWhit^ 12002. for a more thorough discussion of 
cluster mass definitions in dissipationless simulations). It has 
been the common choice by both observers and theorists to 
define the mass within an radii where the average overden- 
sity is greater than a large multiple of a given background 
density, such as Pa(z) and Pm(z)- For low redshift observa- 
tions, the more popular definition has been the Pa iz) as the iso- 
density surface, since no prior knowledge of fi^ is required. 
The question remains what definition is physically more intu- 
itive when comparing across various redshifts. At late times 
(z < 1), clearly the inclusion of the dark energy greatly influ- 
ences the redshift evolution of the critical density compared 
to the mean matter density. For a hypothetical isolated non- 
accreting cluster using the /?a definition will result in the clus- 
ter radius shrinking as time approaches present day. Using the 
RA.m scaling we find that the radial regions at which kinetic 
pressure is in equipartition with thermal pressure and the sharp 
break found in the ICM ellipticity align at ^ 200/?A.m (cf. Fig. 
I2TI ). In BBPS3 we show that the velocity anisotropy has the 
same radial trends as in Figures [T] and QO] and that 200/?A,m 
traces a distinct dynamical region of clusters, the splash-back 
radius, i.e., is caused by the turn-around of earlier collapsed 
shells which minimizes the radial velocity component such that 
the tangential components dominate the velocity. 

D. GAUSSIAN OR LOG-NORMAL SCATTER? 

Previous approaches quantified the scatter around the best-fit 
Y-M scaling relation with a log-normal distribution, i.e. they 
characterized the distribution of SlogYA = log Fa - log Ia. fit 
with a Gaussian. Deviations from this log-normal distribution 
were computed with the Edgeworth expansion, introducing 
substantial higher order moments, such as skewness and kur- 
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Figure 21. The choice for our working definition of virial radius has an i mpa ct on the redshift evolution of both, the kinetic pressure support (left) and ellipticity 
(right) of clusters. The figures shown here are the same as Fig. [T]and Fig. llOl except that the dimensionless radius has been scaled by i?200.m instead of S200- With 
this definition of virial radius, the redshift evolution of both, kinetic pressure support and ellipticity is weakened, especially in the outer regions. 
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Figure 22. Comparison of Gaussian and log-normal scatter relative to the best- 
fit Y-M scaling relation at z = 0. We show the distr ibut ion of the relative linear 
deviation from the mean relation &Y/^/Y (cf. Eq. ll5t with the solid blue line 
and compare it to a Gaussian fit (blue dotted line) and log-normal fit (5 log 7^ , 
red dashed line). The Poisson deviations are shown with the grey band. Here 
we transformed the fit to the <51ogFA distribution into SYi^/Y so they could 
be shown together. The SY^^/Y distribution is fit by a Gaussian better than the 
&\ogY/^ distribution, with X^i^Y/^/Y) ~ 1 and X^ii^ogY/^) ~ 7. Forcing a 
log-normal distribution introduces higher-order moments such as skewness and 
kurtosis as can be seen by the asymmetiic shapes of the tails in the log-normal 
fit. 



tosis (e.g.j Yang et al.ll201 0'). Using a non-linear least squares 
approach we fit a Gaussian to both the (JlogiA and SYj^/Y dis- 
tributions. In Fig. |22] we show that SYj^/Y distribution is a 
better fit by a Gaussian within the (Poisson) uncertainties than 
the (5 log Fa distribution, with ^ 1 compared to ^ 7, re- 
spectively. Hence we suggest to use relative linear deviation 
instead of log-normal scatter for future characterizations of the 
scatter in Y-M relation. 



